Disorder induced transition into a one-dimensional Wigner glass 
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The destruction of quasi-long range crystalline order as a consequence of strong disorder effects is 
shown to accompany the strict localization of all classical plasma modes of one-dimensional Wigner 
crystals at T = 0. We construct a phase diagram that relates the structural phase properties of 
Wigner crystals to a plasmon delocalization transition recently reported. Deep inside the strictly 
localized phase of the strong disorder regime, we observe "glass-like" behavior. However, well into 
the critical phase with a plasmon mobility edge, the system retains its crystalline composition. We 
predict that a transition between the two phases occurs at a critical value of the relative disorder 
strength. This transition has an experimental signature in the AC conductivity as a local maximum 
of the largest spectral amplitude as function of the relative disorder strength. 
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The behavior of interacting electrons in a random en- 
vironment has long been a topic of interest in condensed 
matter physics. When the electrons are confined to a sin- 
gle dimension, many of the complications inherent in the 
study of such a system are vitiated; a number of stud- 
ies have focused on the behavior of the collective modes, 
the various competing manifestations of long-range or- 
der and, of course, the equilibrium structural properties 
of such systems P, S Q • More specifically, in the regime 
in which electronic correlations are expected to be dom- 
inant, the Wigner crystal(WC) phase has been proposed 
as a possible equilibrium configuration in several different 
ID and quasi- ID systemsjl, Q. 

One possible consequence of disorder in any dimension- 
ality is the localization of both electronic wave functions 
and collective modes [1, Q . Recently, the plasmons of a 
disordered ID WC have been reported to exhibit a delo- 
calization transition [Sj]. Subsequent investigations have 
clarified the statistical arrangement of the electrons at 
equilibrium for different types and strengths of disorder 
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In many discussions of ID WC systems at T = a 
primary concern is whether the system is truly a crystal 
in the thermodynamic limit when the effects of quantum- 
mechanical zero-point motion arc included. Previous ef- 
forts have shown that — although there may not be true 
long-range order — the plasmon displacement u{x), cor- 
relation functions decay in space much slower than any 
power law, {[u{x) — u(0)]^) cx y/ln(x). Evidently, this 
implies that strong Bragg peaks would appear in an 
experimental probe of the crystal scattering intensities, 
thus indicating quasi long-range order [1, 0|. Further- 
more, one is justified in the use of classical methods for 
studying the plasma modes, due to the lack of particle 
wave- function overlap at strong unscreened repulsion, ex- 
isting at lower densities. In this letter we focus on the re- 
lationship between the localization of the plasma modes 
in a one-dimensional array of localized charges (a one- 
dimensional Wigner Crystal, or IDWC) and the extent to 
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FIG. 1: Phase diagram of the disordered ID Wigner crystal, 
showing a distinct localized "glass" phase for k < Kc. The 
quantity k quantifies the strength of the Coulomb interactions 
between the localized electrons, in relation to the strength of 
the disorder (see the discussion below Eq. H}). 



which quasi-long-range crystalline order prevails in this 
system in the presence of a "white noise" random po- 
tential at r = 0. We first confirm an vital aspect of 
the plasmon Anderson transition not discussed before, 
in that at and beyond a critical value of the disorder 
strength, measured with respect to the interaction mag- 
nitude, all plasma eigenmodes are localized, in analogy 
to what is seen in the three-dimensional (3D) Anderson 
transition for non-interacting electrons. Strikingly, we 
find that all Bragg peaks at T = also disappear in the 
immediate neighborhood of this critical value of the dis- 
order. This aspect of the transition can be placed in the 
context of earlier work on the Anderson transition, in 
which it is predicted that localization effects tend to sta- 
bilize glassy behavior ll| . The core of these results are 
illustrated by the phase diagram shown in FigUl where 
the phase boundary separates the frequency uj regime 
that the plasma modes of the ID WC system are local- 
ized from the regime in which they are extended. 

The physical system is controlled by the model Hamil- 
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tonian of a standard Jellium Wigner Crystal with an 
added random potential that interacts with each charge 
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with V{x) — A[Y^ a„ cos(27m2;) + &„ sin(27mx)]. The pa- 
ri 

rameters J and A are coupling constants that define the 
dimensionless interaction strength n = J/ A. The random 
variables a„ and 6„ are chosen from different Gaussian 
distributions, yielding a white noise power spectrum with 
a mean of /i = and a variance of cr = 1. The number 
of Fourier amplitudes, N, was set to N = L/A for proper 
scale invariance. We have observed that the critical mo- 
bility edge frequency is generally dependent on the 
ratio N/L of V{x), which represents the average number 
of charges available per potential well. The scaling sig- 
nificance of N in finite sized systems is not immediately 
obvious and must be studied more carefully. 

Since we are interested in quasi long-range crystalline 
order (due to zero-point motion), it is sufficient to study 
large finite systems with periodic boundary conditions. 
The charges were numerically relaxed to their equilibrium 
configuration with the use of methods are outlined in [9] 
and briefly described below. Let us introduce the various 
physical quantities of interest that are needed to interpret 
the results of our investigations. 

First, the central quantity used for studying the 
plasma modes is the dynamical matrix D(R). For 
a finite sized chain of length L, D(R) is an L x L 
symmetric matrix with the structure, D(R— R') = 
Vr'Er" 9mx)/dx^\,^^^_^„ - d^^{x)/dx'\^^^^_^„ 
where </>(x) is defined as the electrostatic potential be- 
tween two charges in a periodic image and the R's are 
the equilibrium positions of the charges, of which are 
not necessarily periodically ordered and must be de- 
termined from a numerical equilibration. The eigen- 
value equation for plasma eigenmodes follows as [12l |: 
TOew2u(R) -I- X^R' D(R- R')w(R') = 0, with w(R) as 
the lattice displacements from equilibrium. Apparently, 
all physical quantities calculated in this article require 
the determination of the equilibrium R's, therefore an 
effective numerical relaxation method is absolutely indis- 
pensable. We have employed a Newton-Raphson proce- 
dure by which one recursively applies the inverse of D(R) 
to the out-of equilibrium particle coordinates R until the 
simultaneous total forces on the charges are sufficiently 
close to zeroQ and consequently the system approaches 
equilibrium. 

An important experimental quantity that follows from 
the plasma oscillations is the AC conductivity. The plas- 
mon propagator G{(jj), is a measure of the response of 
the particle array to an external electric field. G{uj) 
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FIG. 2: (color online) (a) The normalized most extended 
eigenwidth vs. interaction strength showing a crossing point 
at K ~ 0.00013 in (b). (c) Finite size scaling analysis at data 
collapse for l/J « 0.42 ± 0.05. 



can be determined from D(R), by constructing the re- 
solvent: G{uj) =Y.rY.r'V{^'^^- D{R,R')). Applying 
the Kubo formula, the AC conductivity is given by [12l |. 
<j{u!) = iAu)G{ijj). 

As in [8], the lengths are scaled so that the size of 
the region occupied by the charges is unity, (/'(x) = 
7r| csc(7rx)|, and the eigenfunction width is given by 
6 = LV(27r2)^^ ^u,(TO)2u,(„)2sin2(i?„ - i?„), where 
Ui{m) and R„i are the respective ith eigenfunction am- 
plitude and particle position at the site m. 

In terms of the structural properties, we have inves- 
tigated the integrity of the Fourier space of particle 
coordinates, given by: Sk = (X]m=i ^*'^^'")/^' ^ ~ 
21:711/ L, m — 1,2,3... For arbitrarily small disorder 
strength, the eigenfunctions are no longer pure plane 
waves. However, some crystalline aspect to the system is 
conserved, expressed by the standard criterion that there 
exists at least one Bragg peak in the spectrum of Sk- It 
is, in fact, sufficient to track the behavior of the height of 
the first Bragg peak, g{K), as the interaction strength k 
is changed. In the neighborhood of the structural phase 
transition (?(«;) — *■ 0, and the system is best described as 
an amorphous solid or glass. 

Typically in finite sized systems with quenched dis- 
order, various quantities exhibit fluctuations that scale 
like l/VX. For smaller values of L, we can reduce the 
severity of these effects by performing averages over suf- 
ficiently large ensembles. However, at larger system sizes 
these quantities are self-averaging and one can reduce the 
number of required ensembles. Another source of numer- 
ical uncertainty is generated in the relaxation procedure 
outlined in ref. where residual forces can move the 
charges slightly out of equilibrium. Therefore, we have 
relaxed the charges to ~ 10^ in relative force magni- 



3 



tude. 

The essence of the Anderson transition is that aU of 
the eigenfunctions of some Hamihonian, or D(R — R'), 
containing an exphcit form of randomness are locaHzed if 
the normahzed interaction strength, k, is below some def- 
inite value Kc- However, for k > Kc, some eigenmodes are 
localized in a transition that also depends on the value of 
the spectral variable w^. That is, the eigenfunction spec- 
trum is divided by a mobility edge, that separates the 
frequency range in which all the states are localized from 
the set of frequencies for which they are extended d, Q . 
In this particular ID WC system, the plasmon mobil- 
ity edge has been confirmed Q. However, a transition 
involving Kc was not observed. We now report on this 
crucial aspect of the transition, as shown in Fig. [2l Our 
general strategy for determining the existence of such a 
transition is to diagonalize D(R — R'), isolate the local- 
ization length ^ of the most extended of all eigenmodes 
and observe how it behaves as a function of k. We have 
computed ^, normalized by the system size, ^{k,)/L. We 
require f/L = 1 for a truly extended eigenmode. As 
shown in Fig. dja), as k is decreased more states are 
localized until coalesces with the lower band edge, 
signifying complete localization. A bulk transition in the 
thermodynamic limit is confirmed by a common crossing 
point in the un-coUapsed data shown in Fig. [2tb). 

Tests for a transition in the limit of an infinite sys- 
tem have been performed with the use of a finite-size 
scaling analysis, where we assumed the dependence of 
the principal quantities, ^/L = F (^L^^^ (k — Kc)/ Kc) , in 
terms of a universal function F. We have verified that 
the data collapse is controlled by the single exponent 
1/(5 « 0.42 ± 0.05, of which has been optimized by ap- 
plying a analysis in the overlapping regime yielding 
a p value = 0.40, for which the data collapse can be 
interpreted as being in statistical agreement. For the 
specific parameters of our model, we have determined 
Kc « 0.00013 ± 0.00002. We have also expanded earlier 
investigations Q into the strong disorder regime near 
Kc and have observed g{Kc) ~ as shown in FiglJfa). 
Additionally, a focused plot of the critical regime is 
shown in FiglSjb), where g{K) reaches a lower satura- 
tion point precisely at g{nc)- We have verified that this 
is a bulk transition in the limit of a large system size 
as shown in Fig[3]Jd). The disappearance of the first 
Bragg peak is a generic feature of structural phase tran- 
sitions into amorphous, "glasslike" structures. A useful 
visual description is given by FiglH where the relaxed 
particle configurations are plotted along with the ran- 
dom potential V{x), in both the crystal(K — 0.001) and 
amorphous(«; — 0.0001) regimes. Clearly, in the amor- 
phous regime, where n < Kc, the charges tend to clus- 
ter into small crystallite domains inside of the potential 
wells. A characteristic length scale of these domains, Rc 
tends like the effective capacity or width of the individ- 
ual well. It is important that we note that Rc is a disor- 
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FIG. 3: (color online) (a) <?(«:) at different system sizes, (b) A 
close up of g{n) in the amorphous regime where g{n = Kc) —» 
0. (c) L=256, 1st Bragg peak height <?(«;) and most extended 
eigenwidth ^/L, showing a coincidence of the structural and 
delocalization transitions at «: = Kc.{d) g{nc) as a function of 
system size L. 



der dependent quantity and regarding previous work on 
similar systems, Rc can be naively associated with the 
Larkin length [13[ for pinned elastic systems. However, 
it should be noted that we are studying a discrete model 
and by contrast the Larkin length arises as the equilib- 
rium length in a force balancing scheme for a continuum 
model. Therefore, the association of Rc with the Larkin 
length should be taken very loosely, if at all. The re- 
lationship between the plasmon delocalization transition 
and physics of pinning in elastic systems is not known. 

We have confirmed a direct association of the quasi- 
long range crystalline order regime with the existence of 
a plasmon mobility edge by determining a coincidence of 
those respective transitions close to k = Kc- Upon in- 
spection it is also evident from the two most extended 
eigenfunctions of the glassy regime, plotted as ui{R) and 
U2{R) in Fig. (Ha), that the localized plasmon eigen- 
modes tend to be confined near double well potentials. 
Consequently, the total data we have presented develops 
the basis for the phase diagram shown in FigU] 

Let us now focus on the behavior of the cr(w) in the 
regimes separated by Kc- From the definition of (t{uj) 
given earlier, Re[cr(Lu)] cx Im[iujG{uj)] cx p{uj), where 
p{uj) is the plasmon spectral density. It is well known 
from previous studies of Anderson localization that p{uj) 
in the strictly localized regime exhibits, discrete and 
well separated peaks that correspond to the bound state 
spectrumQ. We computed Re[a{uj)] for both k < Kc 
and K > Kc as shown in Fig. |4l Clearly in the local- 
ized regime. Fig. 3] (b) we recover the expected discrete- 
ness, essential to Anderson localization. It follows that 
in the thermodynamic limit, the smaller individual peaks 
would smear out into a continuum. However, this differs 



4 



A 


8 




(b) 




6 






1 \ K=0.001 






K=0.0001 


4 






2 








n 


J 


L 



10 20 „ 40 



0.00006 0.00012 0.00018 
K 



FIG. 4: Re[(j{uj)] (L=128) (a) k > (b) k < Kc, (c) Pre- 
dicted experimental signature of the transition. Largest spec- 
tral amplitude Re[aMAx{io)] forms a local maximum precisely 
at ft/c . 
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FIG. 5: (color online) The distinct structural phases: The 
random potential V(x), plotted together with relaxed equi- 
librium particle coordinates R{x) and the two most extended 
eigenfunctions, ui{R), U2{R), for L=64. (a)K < Kc, in the 
Wigner glass phase, where the particle array has fractured 
into small domains that vibrate with localized plasmon eigen- 
modes. (b)K > Kc, Quasi-long range order with extended 
plasmon eigenmodes. 



tern with short ranged interactions, pinned by disorder 
was studied in where an elastic string was pinned 
into small domains with effective lengths that define the 
Larkin scale Rc- If we consider the AC response of 
the individual domain it has been predicted that the 
peak frequency ujp should scale with the domain size 
ujp ^ 1/Rc l4]. This dependence appears consistent with 
the shift into a spread of larger cOp values as shown in Fig. 
m (b). Intuitively, one would expect a smaller domain to 
vibrate with a higher ujp values. 

We conclude by noting that the results of this paper 
have direct relevance to not only IDWC systems but the 
general properties of bosonic excitations in random me- 
dia and amorphous solids^, The phase diagram 
shown in Fig. [T] is the primary means for linking a sys- 
tems structural properties to the localization of its col- 
lective modes. A more general topic that should be pur- 
sued further is the study of delocalization transitions in 
harmonic oscillator systems in any dimensionality that 
have a non-crystalline configuration at equilibrium. Is 
there a critical power of the interactions between the os- 
cillators that determines whether the model is critical or 
not? For example we have considered a l/R Coulomb 
potential in ID. Thus, for a general power law interac- 
tion I/i?" in any dimension d, what is the critical value 
of a or d for which a Kc should exist? This suggests the 
possibility of multiple parameter scaling features (includ- 
ing N) with a more intricate phase diagram. A broader 
motivation can be extended recent developments in bi- 
ological and optical lattice systems. Another important 
avenue to explore further is the role of the double well 
potential and the various tunneling and energy splitting 
processes. Perhaps the low temperature, universal prop- 
erties associated with the double well potential as gener- 
alized by Anderson-Halperin-Varma theory of amorphous 
solids plays a role in describing the plasmon conductance 
of the glassy regime 17| . Many issues regarding quantum 
mechanical effects such as dissipation and the role of plas- 
mon exchange and tunneling through finite barriers will 
also be addressed in future investigations. 

We thank Professors G. Gruner, S.E. Brown, and S. 
Chakravarty for useful discussions. J.R. Acknowledges 
support of the NSF through grant no. DMR 04-04507. 



from the k > Kc response of Fig. 0] (a), corresponding 
to a large, single and intact domain. Therefore we em- 
phasize that an experimental signature of the plasmon 
Anderson transition would be present in the frequency 
dependent AC conductivity as a shift between these two 
general forms. We have examined the maximum ampli- 
tude which is defined as the maximum value of Re[(j(uj)] 
for a full spectra, as a function of k. Evidently, a local 
maximum forms precisely at shown in Fig. |4] (c). 

Lastly, we mention that cf{uj) for a ID elastic sys- 
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